Global mean first-passage times of random walks on complex networks 
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We present a general framework, applicable to a broad class of random walks on complex networks, 
which provides a rigorous lower bound for the mean first-passage time of a random walker to a target 
site averaged over its starting position, the so-called global mean first-passage time (GMFPT). This 
bound is simply expressed in terms of the equilibrium distribution at the target, and implies a 
minimal scaling of the GMFPT with the network size. We show that this minimal scaling, which 
can be arbitrarily slow for a proper choice of highly connected target, is realized under the simple 
condition that the random walk is transient at the target site, and independently of the small-world, 
scale free or fractal properties of the network. Last, we put forward that the GMFPT to a specific 
target is not a representative property of the network, since the target averaged GMFPT satisfies 
much more restrictive bounds, which forbid any sublinear scaling with the network size. 

PACS numbers: 



Complex networks have appealed a lot of interest in 
the past few years 0, 0, [H, mainly because of the ex- 
tremely broad range of systems that they model, from 
biology to computer science or sociology. Despite this 
variety and their intrinsic topological complexity, many 
real networks have been shown to share some common 
features, such as the small- world property the scale 

free property [J H|, or even fractal scalings Q. A cru- 
cial issue, still under debate, is to understand the impact 
of the topological complexity of such systems on trans- 
port properties. As a paradigm of dynamical processes, 
random walks on complex networks have been intensely 
studied @, H, H, EH ■ an d in particular first-passage 
times have been widely used as a quantitative indicator 
of transport efficiency [HI, [HI- The mean first-passage 
time (MFPT) to a target point was for instance calcu- 
lated in the case of fractal networks [3, [Tp. 

Following the seminal work of Montroll [16| , many pa- 
pers have focused on the MFPT averaged over the start- 
ing point of the walker @, Q3, M, EI W, SH, S3, H, H , 
sometimes called the global mean first-passage time 
(GMFPT) . Recently, a sublinear dependence on the size 
A of the network of the GMFPT to the most connected 
node of a specific network was shown [23|, and was in- 
terpreted as favorable for an efficient trapping. This 
finding, in strong contrast with previously known re- 
sults in the case of regular or fractals [1, [H, [H HH 
lattices, has motivated an increasing number of works 
[H, HI [H, [H, [H that have tried to find examples of 
networks with high trapping efficiency, namely displaying 
weaker and weaker dependence on N of the GMFPT. Re- 
lying on these specific examples, the heterogeneity, and 
more precisely the scale free property was put forward as 
advantageous [H [H| , whereas the fractal property was 
suggested to be unfavorable (24[ ■ 

Here, we propose a general framework, applicable to a 
broad class of networks, which deciphers the dependence 
of the GMFPT on the network size A and provides a 



global understan ding of recent results obtained on spe- 
cific examples [II [Mill, [II [Mill]. We first show on the 
example of a new set of networks that the GMFPT to the 
most connected node can scale as N e , with arbitrarily 
close to despite the fractal property of the network. We 
then present an analytical approach which yields (i) rig- 
orous bounds on the N dependence of the GMFPT, and 
(ii) a simple criterion under which this bound is reached, 
which in particular provides a condition for a sublinear 
scaling with A, which is independent of the scale-free, 
small-world, or fractal nature of the network. Last, we 
show that a sublinear scaling is never representative of 
the network, in the sense that the GMFPT averaged over 
the target site always scales faster than A. 

Definition of the problem and notations. We consider 
a set of graphs {C/ 9 } 9 eN where N g denotes the number of 
sites of the graph Q g at generation g, such that N g — > oo 
when g — > oo. We consider a discrete time random walker 
on Q g . We assume that the transition probabilities Wij 
from site i to site j defining the walk are such that an 
equilibrium distribution P eq satisfying detailed balance 
exists. We further assume that sup X€ g g P eq (X) — > 
when g — > oo. We denote by Fs—,t(T = n) the prob- 
ability that the walker reaches the target site T starting 
from site S for the first time after T = n steps, and write 
T S _ T for the MFPT from S to T. Note that this first 
average 7 is taken over the realizations of the random 
walk. Taking the average of the MFPT over the starting 
point, we define the GMFPT according to : 



GMFPT (T) = <T s ^t)s#t 



S=jtT 

i-^ q (r) 



(1) 



Note that this quantity depends on the target point T. 
Here the space average (•) is taken over the equilibrium 
distribution P^, and slightly differs from the definition 
used in [lj| H2, [H, H3, E3| where the average is taken 



2 



over the flat distribution. It can be checked numerically 
on networks recently studied in the literature that both 
definitions lead to the same scaling with N g . 

Efficient trapping on a fractal network. We first ex- 
hibit a set of fractal networks which extends the so-called 
(u, u)-flowers introduced in (28|, and whose GMFPT to 
the most connected node scales as N e , with 9 arbitrar- 
ily close to 0. The first generation of the graph consists 
in two nodes connected by one link ; then, at each it- 
eration, every link is broken and replaced by k paths of 
u > 2 links ( cf. figure []} . It is clear that this network 
is fractal with a fractal dimension df = ln(fcu)/ ln(u), 
since the diameter of the network at generation g is 
L g ~ u 9 while the number of sites is N g ~ (ku) 9 (the 
usual (u, u) flowers of (HI correspond to the special case 
k = 2). Taking the target as one of the initial node, it 
is easily seen that the determination of the GMFPT is 
actually a ID problem since all the points n(r) at the 
same distance r of the target are equivalent by sym- 
metry, and thus lead all to the same T(r). Noting 
next that for all r e [l,u 9 - 1], P cq (r)n(r) = 2P cq (T) 
and P eq (u 9 )n(u 9 ) = P eq (T), and using the classical ID 
expression T(r) = r(2u 9 — r) (29[, we obtain the fol- 
lowing exact expression: GMFPT(T) = 2 " 9 ( 2 ^ a + 1 ) ^ 
^M^VM***). i n other words, for k large enough, the 
trapping at the hub is arbitrarily efficient on this net- 
work despite its fractal property. 




FIG. 1: A fractal network leading to efficient trapping by the 
hub T (on the most connected sites): case of k — 3, u — 2 at 
generation 5 = 3. 



Lower bound of the GMFPT. In order to gain un- 
derstanding in the real parameters relevant to the scal- 
ing of the GMFPT with the size N, we now derive 
a general lower bound for the GMFPT. This deriva- 
tion follows from the generalization of the Kac for- 
mula [13, HH which we briefly recall here for the sake 
of selfconsistency. We start from the following back- 
ward equation satisfied by Fs->t for n > 2 (see [29j|): 
Fs^rin) = J2j^T w Sj Fj^t(ji — 1), which is completed 
by Fs^rin = 1) = w$t- Laplace transforming and av- 
eraging this equation over S (with a weight P cq (S) as in 



El-©) yields the generalized Kac formula 



P eq (T) 
1 - Pc q (T) 



(s) 



where F 



S->T 



(s) 



n=l 



= (e- s -l)(F 5 _ tT (s))s#T, 
(2) 

i~ sn Fs— >t(h>). This very general 



equation, derived in a similar form in [3l|, relates the 
distribution of the first return time to a site T to the 
distribution of the global first-passage time to T. Ex- 
panding Eq. @ to first order in s yields the classical Kac 

formula T t ^t = l/P cq {T) [HElj]. 
order in s gives : 



In turn, the second 



GMFPT(T) 



1 P cq (T)T 2 T ^ T - 1 

2 1 - P cq (T) ■ 



(3) 



Using next T 2 t^t > ^t^t anc ^ the classical Kac for- 
mula, the above exact expression gives a lower bound for 
the GMFPT: 



GMFPT(T) > 



2P cq (T) ' 



(4) 



Note that this lower bound is in close analogy with the 
one obtained in (32| in the context of continuous space 
Pearson random walks in confinement. 

We now discuss under which conditions this lower 
bound is reached. Strictly speaking, this requires the 
very restrictive condition that the variance of Tt^t is 
zero. More generally we can discuss under which condi- 
tions the right and the left hand side of Eq. ^ share the 
same scaling in the large size limit. To do so, we consider 
a sequence of target sites {T g € <5 g } g eN, which can be for 
instance hubs of the networks at each generation as in 
refs [l^, [22], [23[ ■ Using ^ , and recalling that we have 
assumed P cq (T g ) — > for g — > 00, wc define the minimal 
scaling of the GMFPT for g -> 00 by 



GMFPT (T g ) = 0(l/P eq (T g )) 



(5) 



Eq. ([3]) then shows straightforwardly that this minimal 
scaling is realized as soon as the reduced variance of the 
first return time is finite in the large size limit, namely: 



(T 2 t->t - T T _ T )/T T ^ T = 0(1). 

We now show that this condition for a minimal scaling 
with the network size N g is actually equivalent to the 
transience property of the random walk at the target site 
T g in the large size limit. We first derive an alternative 
exact expression for the GMFPT. Let us introduce the 
pseudo-Green functions [3, HH, [13] defined as: 



H 



(P s ^ T (n) - P cq (T)) . 



(6) 



where Ps—,T(n) is the propagator, namely the probability 
that the walker is at T at time n starting from S. It can 
be shown (see SHU!) that the MFPT is then given by 
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the exact expression: T$->t = P jj* (Ht-*t — Hs^t)- 
Making use of the relation P cq (S)H s ^T = P ca jT)H T ^s, 
which follows from detailed balance (see also Q), we ob- 
tain a second exact expression for the GMFPT: 

This equation provides an alternative condition under 
which the minimal scaling is realized, given by Ht—,t = 
0(1) in the large g limit. From the definition ^ of 
Ht^Ti this condition states that the random walk is 
transient at site T g in the limit g — > oo, i c that in this 
limit, a random walker returns on average only a finite 
number of times to T g [34[ . Conversly, Eq. f7]) indicates 
that if the walk is recurrent at T g for g — > oo, that is 
if HT g ^T g diverges for g — > oo, then GMFPT (T g ) grows 
faster than l/P eq (T g ). 

The lower bound ^ and minimal scaling |(SJ) for the 
GMFPT obtained above call for comments, (i) First, our 
analysis puts forward a very general criterion to have a 
minimal scaling of the GMFPT with the network size, 
namely the type (transient or recurrent) [34[ of the ran- 
dom walk at the target site. We stress that this criterion 
is independent of the scale-free, small world or fractal 
properties of the network. Note that for a generic set of 
graphs {C? 9 } flg N, the type of the random walk for g — > oo 
is a site dependent property [1, [28|, [34| . (ii) Second, the 
minimal scaling §E§ is fully determined by the equilibrium 
distribution at the target site, which is generally much 
easier to obtain than dynamical quantities, and which 
crucially depends on the connectivity of the target site. 
Let us take the classical example of an isotropic random 
walk, for which Wij = 1/fcj if i and j are neighbors and 
else 0, where fc; denotes the connectivity of site i. The 
minimal scaling of the GMFPT to a target T g then reads 
N g {k) /hr g , where (k) is the connectivity averaged over 
all sites, (hi) Note finally that in the case of a recurrent 
random walk at the target the minimal scaling is not re- 
alized, but the scaling of the GMFPT can however be 
sublinear if the growth of the connectivity at the target 
is fast enough. In this case the scaling of the GMFPT de- 
pends on the scaling of Ht ;i ~>t ;i , which generally depends 
both on the network and on the target T g . 

It is noteworthy that our analysis provides a compre- 
hensive view of recent papers highlighting a sublinear de- 
pendence of the GMFPT to a hub on different examples 
of networks, (i) In the example of deterministic scale- 
free graph proposed in , the minimal scaling that we 
predict in Eq.([5]) is indeed realized and the transience 
of the random walk at the target site (as defined above) 
is shown in the limit of large size (since the probability 
to come back at the hub in a finite time is null, as can 
be seen from Eq.(36) from [H^ ] in the large size limit), 
in agreement with our approach, (ii) The authors of 
[22l , f23[ | have studied different examples of small world 
scale-free networks (Apollonian networks [22j and (it, v) 
flowers [H) where the GMFPT to the main hub dis- 
plays a sublinear scaling. In these examples the scaling 



of GMFPT(T S ) is strictly faster than our predicted min- 
imal scaling l/P eq (T g ) (and satisfies the upper bound 
given in ©). Our criterion therefore implies that ran- 
dom walks on such structures are recurrent at the target 
site in the large size limit. 

Bounds on the average GMFPT. As demonstrated pre- 
viously, the GMFPT highly depends on the target site, 
especially in the case of scale-free network where the 
connectivity can be very heterogeneous. Therefore the 
GMFPT to a specific target site cannot be taken as a 
general characteristic of the network. Actually, as we 
proceed to show the GMFPT averaged over the target 
site, defined by (GMFPT) = £ T P oq (T)GMFPT(T), has 
scaling properties with N g which can widely differ from 
the case of a fixed target site studied above. The in- 
equality ^ gives straightforwardly the following lower 
bound for (GMFPT) (see aso (U): (GMFPT) > 
Hence, the averaged GMFPT always scales faster than 
N g , and sublinear scalings discussed above are pointwisc 
properties which are never representative of the network. 
This general inequality sheds some light on the result 
obtained by Bollt and ben Avraham [8| in the case of a 
specific network ((1,2) flowers), where the GMFPT av- 
eraged over a fraction of nodes of the network scales 
sublinearly with N g , while the GMFPT averaged over 
all the nodes is linear. Interestingly, we can also pro- 
pose an upper bound for (GMFPT) following [30| . First 
we define (see also 0, 0) the mean commute time as: 
Tij = T.^j + Tj_>j. The quantity can actually be 
bounded using a very useful electrical analogy. Let us as- 
sign a unitary resistance to each link of the graph. Then 
it can be shown (see (35[) that the following general re- 
lation holds Tij = N g {k)rij, where is the effective 
electrical resistance of the network between sites i and 
j. It is then straightforward to obtain that < dij 
where dij is the distance between i and j. Indeed, dij 
is the resistance of a path of length dij between i and j, 
and any parallel paths can only lower the resistance. We 
therefore finally obtain: 

f < (GMFPT) < (8) 

where (d) is the weighted average over pairs of the point 
to point distance dy . 

Importantly, this shows that the scaling of (GMFPT) 
is much more constrained than the scaling of the GMFPT 
for a fixed target. This is particularly striking in the case 
of small world networks for which (d) ~ In N g : hence in 
case of small- world networks with finite (k), widespread 
in nature this shows that (GMFPT) always scales 
linearly with N g (up to log corrections). Note also that 
these bounds ([§]) are compatible with the linear scaling 
of (GMFPT) with N g reported in the case of Apollonian 
networks [lj| and (1,2) flowers The conditions for 
which the scaling of each of the bounds in ([8]) is realized 
can also be discussed. As for the scaling of the lower 
bound, a sufficient condition for its realization is that for 
any sequence of targets {T g S Gg}g£n, the random walk 
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is transient at T g in the limit g — > oo. Note however 
that this condition is not necessary, and the bound can 
be reached for networks having mixed type properties, as 
in the case of (1,2) flowers already mentioned [8(]. As for 
the scaling of the upper bound, first notice that for any 
tree graph, r%j is exactly the distance dij as discussed 
above using the electrical analogy. We conclude that for 
any tree the scaling of the upper bound is realized. In 
particular we find that (GMFPT) ~ N g \nN g for any 
small world tree (see [36[ for an example). 

Additional comments are in order, (i) First, Eq.© 
provides as a by-product an upper bound for the GMFPT 
itself, leading finally to: 

_L_< G MFPT(r)<™. (9) 

(ii) Second, this upper bound for GMFPT inductively 
gives an upper bound of the trapping time in the case of 
a moving target using the Pascal principle (37| . (iii) Last, 
we underline that in the case of fractal networks, charac- 
terized by a fractal dimension df and a walk dimension 
d w Han explicit scaling of (GMFPT) can be obtained 
(see ja, (3). Indeed, using for instance the asymptotics 
of the MFPT between points separated by a distance r 
[l4| and averaging over r, one gets the following scaling 

( N g if d w < d f 

(GMFPT) - I N g \n(N g ) if d w = d f ; (10) 

{ N^^ Hd W >df 



which depends on the type of the random walk (transient 
if df > d w ). 



Conclusion. We have presented a general framework, 
applicable to a broad class of networks, which provides 
rigorous bounds on the size dependence of the GMFPT 
to a target site. We have shown that the GMFPT has the 
same scaling in the large size limit as this lower bound 
under the condition that the random walk is transient 
at the target site. This shows that the type of the ran- 
dom walk (transient or recurrent) is a crucial criterion 
to determine the scaling of the GMFPT, widely indepen- 
dent of its scale free, small world, or fractal properties. 
This study reconciles recent works on GMFPT for ran- 
dom walks on various network examples. Additionaly, we 
have demonstrated that the scaling of the GMFPT to a 
specific target is not a representative property of the net- 
work, since the target averaged GMFPT satisfies much 
more restrictive bounds, which in particular forbid any 
sublinear scaling with the network size. 



We thank E. Agliari and R. Burioni for useful discus- 
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